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/^-PROTECTED VERTICES IN BINARY SEARCH TREES 



MIKLOS BONA 



I Abstract. We show that for every k, the probabihty that a randomly 

f^ selected vertex of a random binary search tree on n nodes is at distance 

r\l k — 1 from the closest leaf converges to a rational constant c^ as n goes 

. to infinity. 

Oh 

< 

CN 1. Introduction 



1.1. 2-Protected vertices in trees. A 2-protected vertex in a rooted tree 

is a vertex that is not a leaf and is not adjacent to a leaf. In social networks, 

protected vertices may represent participants who have, in the past, invited 

,-^ others to join the network, but have not recently done that. This, and 

C^ other applications led to a recent flurry of interest in studying 2-protected 

H vertices in various kinds of rooted trees. See the articles [3], [1] and [5] for 

'— ' some results. 



J> 1.2. Vertices at level k. We generalize the notion of 2-protected vertices 

lO as follows. In a rooted tree, we say that vertex v is at level k, or is {k — 1)- 

protected if the shortest path from v to any leaf of the tree consists of A; — 1 



\^ edges. In other words, the distance between v and the closest leaf is A: — 1. 



So leaves are at level 1, neighbors of leaves are at level two, and so on. In 
(^ particular, 2-protected vertices are those that are at level 3 or higher. 

C^ In this paper, we will study the numbers of vertices at level k in binary 

^7^ search trees, which are sometimes also called decreasing binary trees, and 

J> which are in one-to-one correspondence with permutations as explained be- 

k> low. 

^ Let p = pip2 ■ ■ -pn he a permutation. The binary search tree of p, which 

we denote by T{p), is defined as follows. The root of T{p) is a vertex labeled 
n, the largest entry of p. If a is the largest entry of p on the left of n, and b is 
the largest entry of p on the right of n, then the root will have two children, 
the left one will be labeled a, and the right one labeled b. If n is the first 
(resp. last) entry of p, then the root will have only one child, and that is a 
left (resp. right) child, and it will necessarily be labeled n — 1 as n — 1 must 
be the largest of all remaining elements. Define the rest of T{p) recursively, 
by taking T{p') and T{p"), where p' and p" are the substrings of p on the 
two sides of n, and affixing them to a and b. 
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Figure 1. The tree T{p) for p = 328794615. 

Note that T{p) is indeed a binary tree, that is, each vertex has 0, 1, or 
2 children. Also note that each child is a left child or a right child of its 
parent, even if that child is an only child. Given T{p), we can easily recover 
p by reading T according to the tree traversal method called in-order. In 
other words, first we read the left subtree of T{p), then the root, and then 
the right subtree of T{p). We read the subtrees according to this very same 
rule. See Figure |1.2| for an illustration. 

Because of this one-to-one correspondence between permutations and bi- 
nary search trees, in our discussion, we will use these two kinds of objects 
interchangeably. 

As a warmup, we try a simple probabilistic approach, which will only be 
successful in the cases of A; = 1 and k = 2. In the case of general k, it will 
provide only a rough lower bound, but that lower bound will be useful in 
the following section. In that section, we use an analytic approach which, 
in theory, provides the exact form of the exponential generating function 
Ak{x) of the total number of vertices at level k in all binary search trees. In 
practice, these generating functions will have a large number of summands. 
However, we will be able to describe them in sufficient precision to find the 
growth rate of their coefficients. 

2. Warm-up: A Probabilistic Approach 

2.1. Two Simple Initial Cases. In this section, we enumerate vertices on 
levels one and two. It turns out that these cases are much simpler than the 
general case, and do not necessitate the general method that we will use in 
later sections. 

In order to alleviate notation, let us agree that for the rest of this paper, 
all permutations are of length n. Let X{p) denote the number of leaves in 
the tree T{p), and let E{X) denote the expectation of X{p) taken over all 
permutations p of length n. 

Proposition 2.1. For all integers n > 2, the equality E{X) = ^^ holds. 

Proof. Let p = piP2 ■ ■ ■ Pn- Let 2 < i < n — 1. Then it is straightforward to 
prove, for instance by induction on n, that the vertex corresponding to pi is 
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a leaf if and only if it is smaller than both of its neighbors, and that event 
has probability 1/3. On the other hand, if i G {l,n}, then pi is a leaf if 
and only if it is smaller than its only neighbor, an event of probability 1/2, 
Therefore, if we denote by Xi{p) the indicator variable of the event that pi 
is a leaf, then by linearity of expectation we get 

n 

1=1 

D 

It is perhaps a little bit surprising that the formula for entries at level two 



is just as simple as the formula proved in Proposition 2.1 Let Y{p) denote 
the number of vertices of p that are at level two. 

Theorem 2.2. Let n > 4. Then the equality E(Y) = ^^"^^^ holds. 

Proof. Let a„^2 be the total number of vertices in all decreasing binary trees 
on n vertices that are at level two. Note that if n > 1, then each leaf must 
have a unique parent, and that parent must always be a vertex at level two. 
However, some vertices at level two are parents of two leaves. We will now 
determine the number dn of such vertices, which will then yield a formula, 

(n + 1)! , 

(1) an,2 = dn 

for On, 2; where n > 4. 

Let Pi be a vertex that is at level two and has two leaves as children. Let 
us assume for now that 3 < i < n — 2 holds. Then pi is larger than both of 
its neighbors, and both of those neighbors pi-i and Pi+i are leaves, so they 
are smaller than both of their neighbors, meaning that pi-i < Pi-2, and 
Pi+i < Pi+2- On the other hand pi must be smaller than both of its second 
neighbors, otherwise its children could not be pi-i and Pi+i. This means 
that if out of the 120 possible permutations of the mentioned five entries, 
only four are possible, since pi must be the middle one in size, its neighbors 
must be the two smallest entries, and its second neighbors must be the two 
largest entries. So if Zi{p) is the indicator variable of the event that pi has 
two leaves as children (in which case pi is necessarily at level two) , then for 
i E [3,n — 2], we get E{Zi) = ^ = ^- If i = 1 or i = n, then pi cannot 
have two children. Finally, ifi = 2ori = n — 1, then an analogous argument 
shows that E{Zi) = ^ = j2' Therefore, since Z = Yl^=2 ^« denotes the 
number of vertices that have two leaves as children (and are therefore at 
level two), then by linearity of expectation we have 

n-l -. Ill 

E(Z) = y E(Zi) = 2- — + (n-4)- — = ^^— ^. 
\ ) Z^ \ ^J 12 ^ ^30 30 

i=2 

Therefore, dn = {n + l)!/30, so formula ([I]) implies that 

(n + 1)! jn + iy. 3 

3 30 10 ^ ^ ' 
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which proves our claim. D 

As a vertex in a rooted tree is called 2-protected if it is not at level 1 or 
2, we can now easily compute the expected number E(Protn) of protected 
vertices in binary search trees of size re. We recover the following result of 
Mark Ward and H. Mahmoud [5]. 

Corollary 2.3. For n> A, the equality 

^,^ X lire -19 
E{Protn) = 



30 

holds. 

2.2. Higher values of k. lik > 2, then finding the total number of vertices 
at level k is significantly more complicated. The main reason for this is that 
if A; > 2, then the unique parent of a vertex at level k — 1 does not have to 
be a vertex at level fe; it can be a vertex at level ^, where 1 < i < k. For 



instance, in the tree T{p) shown in Figure 1.2 vertex 3 is at level two, and 
its parent, vertex 8, is also at level two. 

2.2.1. A simple, hut useful Lemma. Let an,k be the total number of vertices 
at level k in all decreasing binary trees at level k. It is then clear that 
o-n,k+i < o-n,k sincc cach vertex at level k + 1 must have at least one child 
at level k. While finding the exact value of an,k is beyond the scope of this 
introductory section, the following lemma will turn out to be useful for us, 
even if its bound is far from being optimal. 

Lemma 2.4. For each positive integer k, there exists a positive constant jk 
so that if n is large enough, then 

0"n,k ^ 

T >lk- 

re • re! 

In other words, for any fixed A;, the probability that a randomly selected 
vertex of a randomly selected decreasing binary tree of size re is at level k is 
larger than 7^. 



Before we prove lemma 2.4, we need a simple notion. A perfect binary 
tree is a binary tree in which every non-leaf vertex has two children, and 
every leaf is at the same distance from the root. So a perfect binary tree in 
which the root is at level i has 1 + 2 + • • • + 2^~^ = 2^ — 1 vertices. 

We will now compute the expected number of vertices pi that are at level 
k for which the subtree rooted at pi is a perfect binary tree. The expected 
number of such vertices is obviously a lower bound for the expected number 
of vertices at level k. 

Let Qk be the probability that for a randomly selected permutation p of 
length 2—1, the tree T{p) is a perfect binary tree (disregarding the labels). 
It is then clear that Qi = 1, and 

(2) Qk+i = 2k+i_iQk- 
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So Q2 = 1/3, and Qs = 1/63. In particular, Q^ is always a positive real 
number. 

Proposition 2.5. Let p = piP2 • • -Pn be a permutation, and let 2'^~^ + 1 < 
i < n — 2^~^ . (In other words, i is not among the smallest 2^~^ indices or 
the largest 2^^^ indices in p.) Let Pf^ be the probability that the vertex pi of 
T{p) is at height k, and the subtree ofT{p) rooted pi is a perfect binary tree. 
Then the equation 

2 

holds for k > 1. In particular, P^ is a positive real number that does not 
depend on n. 

Proof. The subtree rooted at the vertex pi of T[p) will be a perfect binary 
tree with its root at level k if the following two independent events occur. 

(1) The string pu^fciof 2*^ — 1 consecutive entries of p whose middle entry 
is Pi correspond to a binary search tree that is a perfect binary tree, 
and 

(2) all entries in p\i^k] are less than both entries bracketing p[j^fc], that is, 
both pj_2fc-i and Pj^_2'=-i- 

The first of these events occurs at probability Qfc , and the second one occurs 
at probability TQkTYy^k^ proving our claim. D 

So Pi = 1/3, and P2 = 1/30, as we computed in the proofs of Proposition 
Furthe 

^3 = ^3 



2.1 and Theorem 2.2 Furthermore, 

2 1 1 



63 36 2268 



Now we are in a position to prove Lemma 2.4 



Proof, (of Lemma 2.4) Let Vi{p) be the indicator random variable of the 
event that the subtree of p that is rooted at pi is a perfect binary tree whose 
root is at level k. Then it follows from the definition of Pk that 

E{Vi{p)) = Pk. 

If V{p) denotes the number of vertices of p that are at level k and whose 
subtrees are perfect binary trees, then the linear property of expectation 
yields 



E{V{p)) = {n- 2'')P, 



k, 

k-l 



since we do not allow i to be among the smallest 2 indices or the among 
the largest 2 indices. Therefore, the total number an,k of vertices at level 
k in all decreasing binary trees of size n satisfies 

O-n.k ^ / 1 2 \ P/ 



> ( 1 - - 1 Pfc > ' '^ 



n ■ n\ \ n J 2 

for n > 2 ~^^. This completes the proof, since we can set 7^ = Pkf^- D 
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3. The Analytic Approach 

3.1. A System of Differential Equations. In order to determine the 
exact value of Un^k for A; > 3, we turn to exponential generating functions. 
We recall the well-known fact that the exponential generating function for 
the number of permutations of length n, and equivalently, decreasing binary 
trees on n vertices, is X^„>o '^'^ = 1/(1 ~ ^)- 

For A; > 1, let Ak{x) denote the exponential generating function of the 
numbers of all vertices at level k in all decreasing binary trees of size n. Let 
Bk{x) denote the exponential generating function for such trees in which 
the root is at level k. In both Ak{x)^ and Bk{x), we set the constant term 
to 0. Note that Ak{x) = ^„>i ^^3;", so in particular, the coefficient of x'^ 
in Ak is the expected number of vertices at level A; in a randomly selected 
decreasing binary tree of size n. 

Then the following differential equations hold. 

Lemma 3.1. We have Bi{x) = x, and 

B'kix) = 2Bk-i{x) ■ (,^ - Bi{x) - B2{x) Bu-2{x)\ - Bk-i{xf 

ifk>l. 

Proof. Let T be a binary search tree counted by Bk{x). Let us remove the 
root of T. On the one hand, this yields a structure counted by B'j^{x). On 
the other hand, this yields an ordered pair of binary search trees such that 
one of them has its root at level k — 1, and the other one has its root at 
level A; — 1 or higher. By the Product Formula for exponential generating 
functions (see for instance. Chapter 8 of [2]), such pairs are counted by the 
ffist product on the right-hand side. At the end of the right-hand side, we 
must subtract i?fc_i(a;)^ as ordered pairs in which both trees have their root 
at level A; — 1 are double-counted by the preceding term. D 

Example 3.2. Let k = 2. Then Lemma \3. 1\ yields 

\2 



B'^ix) = 2Bi{x) ■ (r^) - Bi{x) 






1 — X 

Therefore, using the equality i?2(0) = 0, we deduce that 
Lemma 3.3. For k > 1, the linear differential equation 



A'k{x) = j^-Ak{x) + B',{x) 



holds. 
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Proof. Let (T, v) be an ordered pair so that T is a binary search tree on n 
vertices, and t; is a vertex of T that is at level k. Now remove the root of T. 
If the root was v itself, then we get a structure counted by B'f^{x), just as we 



did in the proof of Lemma 3.1 Otherwise, we get an ordered pair {R, S) of 
structures, one of which is a binary search tree, and the other one of which 
is an ordered pair of a binary search tree and a vertex of that tree that is 
at level k. This explains the first summand of the right-hand side by the 
Product Formula. D 

Example 3.4. Setting k = 2, we see that A2{x) is the unique solution of 
the linear differential equation 

2x 



Ai.^) = -, Mx) + 



l-x ' ' 1- 

with initial condition ^2(0) = 0. This yields 



X 



2 



A2ix) 



-gX^ + |x^ — x^ + x^ 



(1-X)2 
4. A CLASS OF FUNCTIONS, AND NEEDED FACTS ABOUT INTEGRATION 

In this section, we define a class of functions that will be useful to describe 
our results. 

4.1. A class of functions. Let PL(x) be the class of functions / : R — )• M 
which are of the form 

(3) /(^)=^«,(i_^)^.ln/-_j , 

where the coefficients Oj are rational numbers, while the exponents bi and Cj 
are non-negative integers. Roughly speaking, PL(x) is the class of functions 

that are " polynomials in 1 — x and In ( j^ 1 " . 

A few facts about PL(x) that are straightforward to prove using integra- 
tion by parts will be useful in the next section. We do not want to break the 
course of the discussion for such technicalities, and therefore we will present 
them in the Appendix. 

4.2. The general form of Ak{x) and Bk{x). Now we are in a position to 
determine the general form of Ak{x) and Bk{x) with sufficient precision to 
deduce the asymptotic number of all entries at level k in all permutations 
of length n. We start with Bk{x). 

Lemma 4.1. For all k > 1, we have 

Bk{x) e PL(x). 



Proof. We prove the statement 
X, and we saw in Example 3.2 



ay induction on k. It is obvious that i?i(x) 
that -Bi(x) = X and B2{x) = 2 In (3^ 



3 

2x — ^. So the statement is true for A; = 1 and k = 2. Now let us assume 
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that the statement of the lemma holds for all positive integers less than k. It 



then follows from Lemma 3.1 that the summands of B'^{x) are all in PL(x), 
except possibly some summands of the form Oj • j^ ■ In ( j— ^ 1 , where Oj 

is a rational number and Cj is a non-negative integer. The integral of each 
such summand is in PL(x) by Fact|7.3[ and integrals of the other summands 



(those that are in PL(x)) are in PL(x) by Proposition 7.1 Therefore, as 



PL(x) is closed under addition, our claim is proved. D 

While the power series Ak{x) are in general not in PL(x), the following 
weaker statement does hold for them. 



Theorem 4.2. For all positive integers k, we have 

Pk{x) 
(l-x) 



(4) A,(x) = /^ + /(x), 



where f{x) £ PL(x), and pk{x) is a polynomial function with rational coef- 
ficients that is not divisible by (1 — x). 



Proof. Lemma 3.3 provides a linear differential equation for Ak{x). Solving 
that equation, we get 

fB',{ x){l-x)Ux C 

[l-xf (l-x) 



(5) Mx)= ''''':rs +n ..2: 



where the integral on the right-hand side is meant with as constant term. 
We saw in the proof of Lemma 4.1 that the summands of B'^{x) are all 



in PL(x), except possibly some summands of the form Oj • j^ • In ( j^ J . 
Therefore, the summands of (1 — x)'^B'Ax) are all in PLfx). Even more 



strongly, each summand of {l — x)B'^{x) is of the form aj(l — x) ' In 



with bi > 1. Therefore, Proposition 7.2 implies that the integral of each 



l-x 



summand is of the form {\ — x)^^^^gi{x)+p(j^^{x), where gi{x) G PL(x), and 
pij\{x) is a polynomial function with rational coefficients. As 6^ + 1 > 2, this 
implies that J B'^,{x){l — x)^ dx = (1 — x)'^g{x) -\-qk{x), where g{x) G PL(x) 
and qk{x) is a polynomial function with rational coefficients, and our claim 
is proved. D 

In other words, though A].{x) contains terms in which (1 — x)^ is in the 
denominator, those terms are simply rational functions; they do not contain 
logarithms. This is important since the coefficients of the power series 

ln(l/(l-x)) 
(l-x)2 

grow faster than those of the terms that occur in ^^(x). (In fact, their 
growth is faster than linear.) 

We are now ready to state and prove the main result of this paper. 



fc-PROTECTED VERTICES IN BINARY SEARCH TREES 9 

Theorem 4.3. Let k > 1, and let an,k be the number of all vertices at level k 
in all binary search trees on n vertices. Then there exists a rational constant 
Ck so that 

lim 7 — — iTT = Ck- 

Proof. Let [x'^]H{x) denote the coefficient of x" in the power series H{x). 
Formula Q shows the general form of Ak{x). The second summand on the 
right-hand side of Q is a function / G PL(x). Each summand of / is of the 

form a.(l - xf^ (in ^)", where, crucially, 6. > and c. > 0, while the a. 

are rational numbers. 

It is proved in Theorem VI. 2. of Analytic Combinatorics [Ij, in particular 
in formula (27) on page 386, that if 6j > and Cj > 0, then 

where the Fj are constants. 

In particular, in each summand of /(x), the coefficient of x" is less than 
iC(ln n)/n for some constant K, and as such, it is negligibly small compared 
to n. Therefore, the contribution of f{x) to [x'^]Ak{x) is negligible, since we 



know from Lemma 2.4 that [x"']^fc(x) > 7a; • n for a positive constant 7^. 

Now we turn to the first summand of formula Q for A]^{x). This sum- 
mand, FJ_^\i is simply a rational function. Its numerator, pk{x) cannot be 
divisible by (1 — x), since that would imply that the coefficients of x" in 
Pfci^J g^Pg g^lY smaller than a constant. (See for instance Theorem IV. 9 in 



[T].) That would be a contradiction since we know from Lemma 2.4 that 



x"]ylfc(x) > 'jkn, and the result of the previous paragraph implies that that 

(l-x)2 



means [x J q_^n2 > 1k^- 



So the rational function /\_\i has a pole of order two at 1. It is now 
routine to prove (see again Theorem IV. 9 in [Ij or the discussion that follows 
here) that 

(7) n/^~c,(n + l) 

(1 — x)"^ 

for some constant Cfc. As we have seen that the contribution of /(x) in (E| 
is insignificant, ^ and ([T]) together imply that 

[x^]Ak{x) ~ [x"]-f^ ~ ck{n + 1). 
(1 — x)^ 

Finally, we prove that Ck is rational. In order to see this, note that if n is 
large enough then 



r^ni Pkjx) _ ax + b 

^ J(l-X)2 "^^ ^(1-X)2^ 
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where ax + b is the remainder obtained when pk{x) is divided by (1 — j;)^. 
As pk has rational coefficients, both a and b are rational numbers. However, 

. „, ax + b . „. / a + b a \ , ,,, ^, 

[^ ] (13^ = [^ H (13^ - (13^ J = (" + ')(^ + ^) - «' 

so Cfc = a + 6, which is a rational number. D 



5. Examples 
5.1. The case of /c = 3. Determining the value of C3 requires finding B'^{x) 



first. We can do that by using Lemma 3.1, since Bi{x) and B2{x) have 
already been computed in Lemma |3.1| and Example |3.2[ A routine compu- 
tation that we carried out using Maple leads to 



ln(l/(l-x)) , ,^,^,,,,, ^,, 2 x3 2^ 



(8) B'^{x) = 4^ ^+4xln(l/(l-x)) 



,4 



-X 



2,1-x 3 
l-x ^" ' 3 -^/(l-^))- 9 



(9) - A-^ 41n(l/(l-x))2 + ^ln(l/(l-x)) 



Now we can solve the differential equation provided by Lemma 3.3 with 
fe = 3, to get 

, , , 1721 x^ x^ bx^ 2x^, ,, ,,, ,, 23x^ 

A?.ix) = -, ^TTv \ \ ln(l/(l-x))H 

^ ' 8100(1 -x)2) 81 324 54 9 ^ '^ " 324 

4x3 349x3 14x2 979x2 

- ln(l/(l-x))H \ ln(l/(l-x))H 

45 ^ '^ " 2025 15 ^ '^ " 2700 

82;, ,, .,, ,, 4219x 4x , ,^ ,,^ ,,2 

- — ln(l/(l-x))H ln(l/(l-x))2 

5 ^ '^ " 4050 3 ^ '^ " 

4 . ,. .^9 1721 22 , ,, 
+ - n 1/ 1 - X 2 n 1/ 1 - X . 



It is now clear, by the proof of Theorem 4.3 that 



x"bl3 X 1721 

C3 = lim \ ^ ^\/ = 0.2124691358. 

n + 1)! 8100 



n—^cx3 



5.2. The case of A; = 4. Determining the value of C4 is conceptually the 
same as determining C3. However, the computation becomes much more 



cumbersome. Lemma 3.1 provides a formula for B'^{x) as a sum. According 
to Maple, that sum has 52 summands of the form OjX ' (ln(l/(l — x)))'^'. 



Using that expression for B'^{x), we can compute A4{x) using Lemma 3.3 
Maple obtains a solution that has 59 summands. However, only 17 of these 
59 summands contribute to P4(x), and therefore, only these 17 summands 
influence C4. The value we obtain for C4 is 

250488312501647783 , „ 

C4 = = 0.1091543117. 

2294809143026400000 
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6. Further Directions 

The data that we computed, c\ = 1/3, C2 = 0.3, cs = 0.212, and ca = 
0.109 suggest that the sequence ci, C2, • • • is log-concave. Is that indeed the 
case, and if so, is there a combinatorial proof? We point out that it is not 
true that for any fixed n, the sequence fln.i, an,2, • • • j On,n is log-concave. For 
instance, n = 4 provides a counterexample. However, it can still be the case 
that for every k, there exists a threshold N(k) so that if n > N(k), then the 
sequence o„^i, an,2, • • • , an,fc is log-concave. 

7. Appendix:Needed facts about integration 

Proposition 7.1. The class PL(x) is closed under integration with respect 
to X. 

Proof. We need to show that J(l — x) In f j^ 1 dx G PL(a;). We prove 

this by induction on c, the inital case of c = being obvious. Integration 
by parts yields 

(10)/(l-.)'ln(^)^. = -ln(^)^a^ 

By the induction hypothesis, the integral on the right-hand side is in PL(x), 
proving our claim. D 

A special case of the previous proposition will be particularly useful for 
us. 

Proposition 7.2. Let 6 > and c> be integers. Then 

/(l - x)''ln (y^) dx = {l- x)^+i • g{x) +p{x), 

where p is a polynomial function, and g{x) G P'L(x). The integral on the 
left-hand side is taken with constant term 0. 

Proof. Induction on c, the initial case being that of c = 0. If 6 = 0, then the 
statement is true, since J 1 dx = x = {1 — x) ■ {—1) + 1. If 6 > 0, then the 
statement is true, since /(I — x)^ dx = {1 — x)^+^ • ^^. 

The induction step directly follows from (10) and from the induction 
hypothesis. D 

We will also need the following. 

Fact 7.3. For all non-negative integers c, the equality 

^ -Inf-^l dx = — *— ln(-^) +C 



1 — X VI — 2;/ c+1 Vl 

holds. 
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